(N 



- 

>< 



Expanded Vandermonde powers and sum rules for the 
two-dimensional one-component plasma 

Gabriel Tellez* and Peter J. Forrester^ 



O _ 

, *Departamento de Fisica, Universidad de Los Andes, Bogota, Colombia 
Q_|] ^Department of Mathematics and Statistics, The University of Melbourne, Victoria 3010, Australia 

< 
(N 



Abstract 



The two-dimensional one-component plasma (2dOCP) is a system of N mobile particles 
of the same charge g on a surface with a neutralising background. The Boltzmann factor 
of the 2dOCP at temperature T can be expressed as a Vandermonde determinant to the 
power r = q^/{kBT). Recent advances in the theory of symmetric and anti-symmetric 
Jack polymonials provide an efficient way to expand this power of the Vandermonde in 
their monomial basis, allowing the computation of several thermodynamic and structural 
properties of the 2dOCP for N values up to 14 and F equal to 4, 6 and 8. In this work, we 
explore two applications of this formalism to study the moments of the pair correlation 
^ function of the 2dOCP on a sphere, and the distribution of radial linear statistics of the 

^ , 2dOCP in the plane. 

o , 

^ ■ Key words: Coulomb gas; one-component plasma; Jack polynomials; sum rules. 

^ 1 Introduction 



Generally the two-dimensional one-component plasma (2dOCP) refers to a model in classical 
statistical mechanics consisting of mobile point particles of the same charge q and a smeared 
5-H out neutralising background, confined in a two-dimensional surface. All the charges (point and 
- -' continuous) interact via the solution v{f,f') of the Poisson equation on that surface. For the 
plane this is 

v{r,r ') = —log (^\r — r'\/ , (1.1) 

where L is an arbitrary length scale, which will be set equal to one from now on, while for a 
sphere of radius R, 

vi{9, 0), (9', 00) = - logi2R\u'v - W|), (1.2) 

where with (^,0) and (6'',0') spherical coordinates, M,f (and similarly m', f') are the Cayley- 
Klein parameters 

M := cos(0/2)e^'^/2^ w := -isin(^/2)e-^"^/2 (13) 
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(see e.g. |9i §15.6.1]). In the latter case, due to the rotational invariance of the sphere, the 
background-particle interaction only contributes a constant, and the Boltzmann factor, at tem- 
perature T, is then computed to be 



\2RJ 



i<j<k<N 



where F = g^/ {ksT). In the former case, the particles couple to the background via a harmonic 
potential towards the origin, and the Boltzmann factor is then (see e.g. [9l eq. (1.72)]) 

^^e--rp,Ef.,k^.P/2 -Q 1^,-^,1^, = ^-TN^i{i/2)Xo^R-m ^ (1.5) 

where ph is the density of the background. Although the derivation of (11. 5p requires that the 
mobile particles and background be confined to a disk of radius i?, for our purposes below 
we will relax this constraint by allowing the particles to locate anywhere in the plane. This 
situation will be called the soft disk geometry. 

Both (II. 5p and (ll.4p have significance in a field of theoretical physics quite distinct from that 
of classical plasmas. Consider for definiteness (II. 5p . According to the well known Vandermonde 
determinant identity (see e.g. [9l eq. (1.173)]), for F = 2 and after introducing polar coordinates 
f = (r, 6) and then complex coordinates Zj = rje^^^ , this can be written in the Slater determinant 
form 



A, 



det[0,_i(^fc)],,fe=i,...,jv] , ^j-i{z) = z^-^-^P'"-^ (1.6) 



Moreover, the single particle wave functions (j)j-i{z)^ j = 1,...,N are, with rrpb = 1/4Z^ 
/ = ^JhcjeB the magnetic length, ground state (lowest Landau level) eigenfunctions of the 
Schrodinger equation for particles confined to the plane in a perpendicular magnetic field of 
strength B (see e.g. [HI §15.2.1]). The integer j has the interpretation of being proportional 
to the centre of the corresponding classical cyclotron orbits. Thus (11.61) is the absolute value 
squared of iV non-interacting fermions in the plane, subject to a perpendicular magnetic field 
which in turn is required to be sufficiently strong that the spin degree of freedom of the fermions 
is frozen out. For F/2 an odd integer, it is a celebrated result of Laughlin [15], that the (un- 
normalized) wave function 



(det[0,_i(.,)],,,=i,...,^])'^' = e-'^''^^S.-l^^^l^/2 n (^^-%)'^' 



is an accurate trial wave function for the fractional quantum Hall effect in the case of filling 
fraction 2/F. Furthermore, the same holds for 



n ^, 

^<j<k<N 



in relation to fractional quantum Hall states on the sphere |18j . 

We have highlighted the plasma system in the plane, and on a sphere. As is usual in 
statistical mechanics, bulk properties such as the free energy per particle are expected be 
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independent of the boundary conditions. In two-dimensions, the sphere is distinguished by the 
usual correction to the bulk free energy per particle — the surface tension — being absent due 
to the homogeneity of the geometry. This feature of the sphere makes it well suited for studies 
in which the aim is to extrapolate bulk values from finite data [HI [IDl [IS]- Extrapolation 
of finite N data for the 2dOCP is to be a primary concern of this present paper, and as such 
we too will make use of spherical geometry and consider the plasma system specified by the 
Boltzmann factor (11. 4p . 

Specifically, we aim to probe the finite analogues of the known sum rules for the bulk 
two-particle correlation function P{2)(^, 0), or more conveniently the truncated bulk two-particle 
correlation function p^)(r, 0) := P{2){r, 0) — p^, where p denotes the particle density. These sum 
rules read [TTl | 



pl){f,6)df=-l (1.7) 
rVf2)(r,0)(ir = -4 (1-8) 



p/ rV;„(r,0)dr=-^(^l--) (1.9) 



It is straightforward to show from the definitions that (11.71) is exact for all on the sphere. In 
fact the sphere plays no specific role; replacing p on the LHS by p(0) this holds for a general one- 
component system. Less obvious and special to the sphere geometry is that (11.81) . appropriately 
modified, can also be made exact for finite N. The derivation is given in Section [3] below. 

To probe (11.91) and (ll.lOp for finite we have available an exact expression for p^)(^, r') at 
the special coupling F = 2 [5j. More generally, for F an even integer (11.41) is the absolute value 
squared of a multivariable polynomial. In the case that F/2 is also even, a recent advance [3] 
identifies the polynomial in terms of the symmetric Jack polynomials (see [H] and ^ for text- 
book treatments), while for F/2 odd an even more recent paper [1] shows that the polynomial is 
an anti-symmetric Jack polynomial. These Jack polynomials have certain structural properties 
(related to their monomial basis expansion) which makes them the most efficient known way to 
carry out exact numerical computations at these couplings. We review the formalism of such 
computation methods in Section [21 and we give too an alternative derivation to the result of [4j, 
[21] for the eigenoperator of anti-symmetric Jack polynomials. The results of calculations based 
on this formalism are reported in Section [3l In Section HI an approximation to the moments is 
probed using a formalism based on the direct correlation function. 

One viewpoint on the moments of p^^ (r, 0) is that they are the averaged value of the linear 

statistic Ylh=i the system perturbed by a particle being fixed at the origin, minus the 

average value of this same linear statistic before the perturbation. This suggests the question 
of computing fluctuation formulas for linear statistics, and we take up this task in Section [5l 
Section [5] is independent from Sections [3] and HI Some concluding remarks are made in Section 

m 
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2 Preliminary material 



In this section the expansion formulas underlying the exact calculations of the moments — 
both analytic and numeric — will be presented. Before doing so, for the sphere geometry, it is 
convenient to map the particles to the plane by applying a stereographic projection. With the 
latter carried out by mapping from the south pole to the plane tangent to the north pole, and 
with {9, 0) the usual spherical coordinates, this is specified by the equation 

g 

z = 2Re"^ tan-, z = x + iy, 

and we then have 



N 



^<j<k<N 
N 



^<j<k<N 



As previously remarked, for F even the product of differences in both (11. 5p and (11.41) is 
the absolute value squared of a polynomial. This remains true of the stereographic projection 
of (II. 4p as seen in (12. ip . In the case T = 4p the polynomial is symmetric, while in the case 
r = 4p + 2 the polynomial is anti-symmetric. As a consequence, the two case need to be treated 
separately. 

2.1 The case T = 4j9 

Following [19], let /i = (/ii, . . . , /xat) be a partition of pN{N — 1) such that 

2p{N - 1) > /ii > /i2 > ■ ■ ■ > /ijv > 0, (2.2) 

and, with rrii denoting the corresponding frequency of the integer i in the partition, define the 
corresponding monomial symmetric function by 

m^{z^, . . . , zjv) = YTzn E <(i) ■ ■ ■ <iNy 

We expand 

n (^;.-^.)'^ = E^r(2pK(^i'---'^'^)- (2-3) 

The significance of knowledge of {c^^''(2p)} in this expansion is that a simple calculation shows 
m eq. (2.7)] 

= iV!.-E(£pM!nG„. (2.4) 

a lli"^*- 1=1 
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where 



/•oo 

GM--='^ / r'^^''9{r^)dr. (2.5) 



2.2 The case T = 4j9 + 2 

We again follow [19]. We now take to be a partition of {p + 1)A^(A^ — 1) such that 

(2p+l)(Ar-l)>/ii>^2>--->/U^>0 (2.6) 

(note that the strict inequalities between the parts of the partition implies = 1, in the 
notation used below (12. 2p ). We then expand 

n {z,-z,f'^^' = Y,cf\2p+l)A{z^,^...z%-). (2.7) 

where A denotes anti-symmetrization. From the definition of the Schur polynomials s^{zi, . . . ,Zn 
and with ^at := (A^ — 1, — 2, . . . , 0), this is equivalent to the expansion 

n {z,-z,f^ = J2i''\2p+l)s,^sA^,,...,zr,), (2.8) 
and now for the integral (12.41) we have 



N 



/;v,rM = iV!7r^$^(c!f)(2p + l))^n^Mr (2-9) 



=1 



2.3 The coefficients c\f\2p) and c\f\2p+l) 

It has recently been observed [31 H] that the products in (12. 3p and (12.81) can be expressed in 
terms of the Jack symmetric polynomial P!f\z) and Jack anti-symmetric polynomial Slt\z) 
respectively, where z := {zi, . . . , z^). Generally the nonsymmetric Jack polynomials are eigen- 
functions of the differential operator [9l eq. (11.63)] 



j=l "'^ j=l ^ 

where interchanges Zi and zj. The operators of symmetrization and anti-symmetrization 
commute with this operator, and so with Mjk = 1 and Mjk = —I respectively, (I2.10p is the dif- 
ferential eigenoperator for Pk°'\z) and Slt\z). This characterisation specifies the polynomials 
uniquely when supplemented by the structural formulas 

Pi")(z) = m,(-2) + ^c,p(a)m,(z) (2.11) 
Sl^\z) = s..sA^)+ c.,(«)s,_5,(z), (2.12) 

p— <5]v<K— <5]v 



where in f l2.12p all parts of /t are required to be distinct, and fi < k refers to the dominance 
partial ordering on partitions < specified by Yli=i — Yli=i '^i^ {k = 1, . . . , N). 

We remark that in the work |3] the anti-symmetric polynomials, defined as in [T] as the anti- 
symmetric Jack polynomials, where not considered directly. Rather attention was focussed on 
Y[i<i<j<Ni^i~^j)'^'^"\^)' which was shown to be the polynomial eigenfunction of the differential 
operator 

Simple manipulation reduces this to 



Zi Zj ^ 



ZiZj 



J=l J 

Comparison with fl2.1Up reveals that up to a constant this is equal to H^^'^^^ with a replaced 
by a/{l — a) and Mij = —1, telling us that fl2.13p is an eigenoperator for Slt^^^~°'\z). This 
conclusion is consistent with the known relation jl] 

l<i<j<N 

The products in (12. 3 p and (12. 8 p relate to the Jack polynomials with a negative parameter. 
Explicitly we have pj (see ^ for a derivation based on earlier literature) 

n iz,-z,)'^ = P,,sA^;-2/i2p-l)) (2.14) 

^<j<k<N 

and [1] 

H {zk-z,yP+' = P^,,+,)sA^;-2/{2p+l)), (2.15) 

i<j<k<N 

where generally for c G Z"*", cn := (cki,ck2, • • • ,ckn)- The computational advantage of these 
formulas in computing the coefficients c^^"* (12. 3p and (12. 8p is that the partitions /i for which 
c}/^'* are nonzero is greatly restricted by the dominance ordering exhibited in (12. lip and (I2.12p . 
Furthermore the eigenfunction characterization allows recursive formulas for the coefficients in 
these formulas, and thus the coefficients {c|j^^} to be obtained. 

In relation to (12. lip , let p = (pi, p2, • • • , Pn), and let p be constructed from p by first adding 
r to Pi and subtracting r from pj with the requirement that upon suitable reordering to be a 
partition it satisfies p < p < k. Also define 

^ 2 
€«(«) := > ^ - 1 («-!))• 

i=l 



6 



Then ^ §VI.4 Example 3(d)] 



^-^= e.(a)-e» ! E ((/^^ + ^) - (z^^' " ^))^-(«)- (2-16) 



Regarding now ( I2.12p . let p and /x be as above and define 

^ 1 
el{a) := 5^ + 2^(1--)). 



Further define (— l)^sw the sign of the permutation needed to be applied to p so that it is 
a partition. Then we have [H [21] 



3 Moments of the correlation function on the sphere 
3.1 Second moment 

The partition function of the 2dOCP on a sphere is [51 [12] O §15.6] 

Q = j^J e'^''dS^...dSN (3.1) 

n i^fc-^rn^ 

l<j<k<N j=l ( 1 + 



iV!(2/?)^^r/2 y^,, l^'^ ^^-1" n ^ i^^i, ^'2+(iv-i)5 ' (3.2) 



where U is the potential energy of the 2dOCP and zj = TjC^"^^ are the coordinates of the 

stereographic projection of a sphere of radius R. Equivalently, one can factor out the radius, 
introducing Zj = Zj/{2R), to obtain 

Q = ^/ TT \~Zk-z/Y\ ^ F- (3.3) 

iV!(2/?)^^(l-) A- JXn M (1 + |%P)^+^^-^)^ 

Let us set 

AT r 

Qi = ^ e-^^(^i'-'^— (3.4) 



where the particle number is fixed at the north pole. By definition, the density at the north 
pole is 

P(i)(0) = ^ = Pf>- (3.5) 

In the coordinates of the stereographic projection of the sphere of radius R, Qi can be written 

as 



l<j<k<N-l j=l I _|_ 



(2i?)2 



or, using the scaled coordinates 

Ql = 7 / TT 1^^ - %r TT ' / .r (3.7) 



Computing the derivative of Qi with respect to R, using expression (13. 6 p gives 
R dQi N^V 



Q dR 2 
Here the last integral appears from the derivative 



P(i)(0) + 2 (^2 + ^(iV - 1)^ j p(2)(e) sin^ ^ dS. 



(3.8) 



2+(Af-l)£ 

d I I \ , , |;z,f/(2i?)^ , , 

TT^ ^ 2(2 + (iV - 1)^) ' : 3.U^ (3.9) 



(2R) 
2+(7V-l)§ 

r, / 1 \ . 



2(2 + (iV-l)-) --^ sin^^. (3.10) 



On the other hand, starting with expression (13. 7p gives 



^^% = 2({N-1)-N^)p^,^{0). (3.11) 



Q dR V 4 

Equating (13. 8p and (13. lip , using the fact that P(i)(0) = pb, and subtracting on both sides 
/ plsm'^{e/2) dS = Npb/2, gives 

where p^)(^) is the truncated pair correlation function between the north pole and a point 
located on the sphere at {9, 0). 

In the planar limit R — )■ oo, N — )■ oo, with pb = N/IAttR"^) finite, the length of the chord 
r = 2_Rsin(6'/2) becomes the distance between the two points, and one recovers the Stillinger- 
Lovett second moment sum rule (11.80 . 

3.2 Higher order moments 

Let us define 

l2n = Pt j {2Rsm{e/2)f-h{e)dS, (3.13) 



where h is the total correlation function 



m = p[2){Q)ipi- (3.14) 
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Using the formalism presented in Section |2] and [19], the density at a given point on the sphere 
is 

N ^ N 

rti,! I /V — M-^ -I- 11! . — . I r): 'I I /vn^ 

P(1)W 



X 



where 6 is the angle from the north pole of the sphere, x = tan(^/2), and 



^sphere 



TV 



r 



-/i,!((iV-l)f-^,)! 

(3.15) 



(3.16) 



is (up to a multiplicative constant) the partition function of the 2dOCP on the sphere. Since 
the sphere is homogeneous, the density does not depend on x and should be equal to pb. This 
gives an interesting relation satisfied by the coefficients, (c|/^''(r/2))^ 



N{1 



■^sphere ^ 1 Ij I'H- 



N 



((^-i)o-/^0!E 



X 



^-/ifc!((7V-l)^ -/.,)! 

(3.17) 

Putting a; = 0, then relates a sum of all admissible partitions and those with = 0, 



^E^^^Vr^n«'((iv-4-«)!=((iv-i)i+i)! E ^Til^n«K(iv-i)i-.,)! 



The two-point correlation function is given by 



with llat—O 



(3.18) 



N-l 



N-l 



with /^^— 



From these expressions, we obtain the 2n-moment of the total correlation function 



(3.19) 



'2n 



N{{N -1)^ + 1 + n)\Z,- 



X 



^sphere 

V- (cr(r/2))^ /^^!((iV - l)f - /i/)! y^^ (/X. + r.) 



AT 



with fij^—O 

In the case n = 1, the sum 

N-l 



k=l 



n + 1 



.(3.20) 



E 

fc=i 



(/ifc + n)! 



^-1 /T/v \ 

5^(/i. + l) = (iV-l) — + 1 



(3.21) 



simplifies and it is independent of the partition fi, which lead to the sum rule discussed in the 
previous section. 



r - A^r - 4 
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(3.22) 



Unfortunately, for n > 2, no such simplification seems possible. 

For r = 2, only one partition appears in the expansion and the result for any n is 



'2n 



-AT- 



n\N\ 



2. 



{N + n)\ ' 

In the thermodynamic limit, — )• oo, the values of I 4, and Iq are also known for any F 



(3.23) 



h 
h 



F-4, 

^(F-6)(8-3F). 



(3.24) 
(3.25) 



For finite A^, we computed the moments numerically from (13.201) . Since the maximum value 
of N for which the numerical calculations where possible is small, the results for the moments 
show important finite size corrections with respect to the known values for — t- 00. The exact 
results for n = 1 (any F), and the one for F = 2, suggest that the finite size corrections can be 
understood as a series expansion in powers of Indeed, for n = 1, 



k=i ^ ^ 



(3.26) 



and for F = 2, 



n{n + 1) _^ n{n + l){n + 2)(1 + 3n) 



2A^ ■ 24A^2 
n^{n + lf{n + 2)(n + 3) n{n + l)(n + 2){n + 3)(n + A){lhn^ + 30n2 + 5n - 2) 



48iV3 



5760iV4 



+0{l/N^ 



(3.27) 



Therefore, we decided to fit the numerical data obtained for the moments with a expansion. 
The results are shown for F = 4 in tables 13.11 13.21 13.31 for F = 6 in tables 13. 4[ 13. 5[ 13. 6[ and 
for F = 8 in tables 13. 7[ 13.81 13.91 For F = 4, the coefficient of order (which gives the value of 
l2n when — )• 00), shows an acceptable convergence to the expected value for and /e, see 

m 



tables 13. H 13.21 Also, from table 13.31 (Jg, F = 4), the zero order coefficient seems to converge 
to a value close to —30. This allows us to obtain an estimate for the value of 8-th moment, 
for F = 4, in the thermodynamic limit, Jg ~ —30. Unfortunately, for F = 6 and F = 8, since 
the maximum value of A^ for which the calculation were possible was smaller and the finite size 
corrections still very large, we were not able to obtain any estimate for Jg. 
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N 


h 


a + 


b/N + 


c/N^ + 








a 


b 


c 


d 


2 


-1.06666666666667 










3 


-0.73469387755102 










4 


-0.552915766738661 










5 


-0.437781621713968 


0.076709 


-2.81362 


1.30721 


-0.506965 


6 


-0.361584090880502 


-0.042858 


-1.37881 


-4.31243 


6.66705 


7 


-0.307439760694233 


0.0121923 


-2.20457 


-0.238709 


0.0610148 


8 


-0.267158562552772 


-0.00224307 


-1.94473 


-1.7833 


3.09245 


9 


-0.236094785664912 


-0.00181613 


-1.9537 


-1.72096 


2.94899 


10 


-0.211435346122364 


0.0000393273 


-1.99823 


-1.36657 


2.01384 


11 


-0.191399568920847 


0.000150015 


-1.99312 


-1.41239 


2.15017 


12 


-0.17480728582518 


-0.000241804 


-1.99036 


-1.43984 


2.24104 


13 


-0.160846001528985 


-0.000138037 


-1.99379 


-1.40227 


2.10407 


14 


-0.148938895568825 


-0.0000794746 


-1.99589 


-1.37703 


2.00357 














oo 














Table 3.1: 1/N expansion coefficients for when F = 4. 



N 


76= a + b/N + c/N^ + d/N^ 






a 


b 


c 


d 


4 


5.1605471562275 










5 


6.11292630320115 










6 


6.56602660174796 










7 


6.78048920359037 


5.49915 


27.6653 


-150.595 


138.066 


8 


6.87498575554362 


5.58979 


26.0337 


-140.897 


119.031 


9 


6.9075123159669 


5.74333 


22.8094 


-118.48 


67.4414 


10 


6.90772192712509 


5.87714 


19.5982 


-92.9236 


0.00502236 


11 


6.89103601184945 


5.92331 


18.3514 


-81.7487 


-33.2425 


12 


6.86574608931141 


5.94863 


17.592 


-74.1805 


-58.301 


13 


6.83613198600859 


5.96717 


16.9702 


-67.3591 


-83.1735 


14 


6.80567000027031 


5.97905 


16.5534 


-62.3686 


-103.044 














oo 
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Table 3.2: 1/N expansion coefficients for Iq when F = 4. 
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N 


/g = a + b/N + c/N^ + d/iV^ 






a 


b 


c 


d 


6 


40.3968590167648 










7 


33.7964087943161 










8 


26.9657510022336 










9 


20.6758161117942 


-49.3234 


709.717 


-95.7946 


-5595.45 


10 


15.1585480309705 


-38.7061 


454.9 


1932.12 


-10946.6 


11 


10.414558865846 


-33.7202 


320.28 


3138.71 


-14536.5 


12 


6.3618885422448 


-31.5027 


253.756 


3801.74 


-16731.8 


13 


2.89953248501109 


-30.4953 


220.513 


4166.41 


-18061.5 


14 


-0.0685490396317987 


-30.1108 


206.671 


4332.12 


-18721.3 



Table 3.3: 1/N expansion coefficients for /g when F = 4. 



N 


h = a + b/N + c/N^ + d/N^ 






a 


b 


c 


d 


4 


1.77112299465241 










5 


1.92727455514225 










6 


1.97167595506035 










7 


1.99917180579183 


2.90767 


-14.9074 


84.3562 


-171.645 


8 


2.01319334170647 


1.72398 


6.39891 


-42.2979 


76.9278 


9 


2.02020128573558 


1.95459 


1.55616 


-8.62918 


-0.556218 


10 


2.02414460391682 


2.08557 


-1.5874 


16.3883 


-66.5709 


11 


2.02625080279343 


1.98854 


1.03254 


-7.09409 


3.29402 


12 


2.02719976576841 


1.98147 


1.2446 


-9.2077 


10.2923 














oo 


2 











Table 3.4: 1/N expansion coefficients for when F = 6. 
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AT 


76= a + b/N + c/N'^ + d/N^ 






a 


7 

b 


c 


7 


A 

4 


29.2857260386672 










r 




24.3553395496129 













19.3570221116088 










1 


15.8157978386411 


23.4057 


-367.801 


3052.6 


-5949.24 


o 

o 


13.2137920826172 


-1.67692 


83.6859 


368.754 


-681.887 


y 


11.2455089383337 


-0.844904 


66.2135 


490.229 


-961.446 


iU 


9.75173060605997 


4.89476 


-71.5384 


1586.5 


-3854.23 


11 


8.58562698480205 


0.124033 


57.2712 


431.99 


-419.314 


12 


7.65247243254821 


-0.286731 


69.5941 


309.171 


-12.6574 














OO 













Table 3.5: 1/N expansion coefficients for le when F = 6. 


N 


/g = a + b/N + c/N^ + d/N^ 






a 


6 


c 


d 


4 


157.436163836164 










5 


-33.5270610707337 










6 


-196.111013214468 










7 


-293.318966661244 










8 


-343.529970425059 










9 


-366.832915368496 


129.267 


-11229.9 


71322.9 


-93943.9 


10 


-374.243763793112 


161.032 


-11992.2 


77390 


-109954 


11 


-373.310338904328 


-86.8351 


-5299.81 


17406.2 


68510.9 


12 


-368.298029500803 


-183.485 


-2400.33 


-11492.1 


64194 


Table 3.6: 1/N expansion coefficients for /§ when F = 6. 


N 


/4 = a + b/N + c/N^ + d/N^ 






a 


b 


c 


d 


3 


5.9504132231405 










4 


5.35216744227873 










5 


5.31949584221798 










6 


5.03608123916663 


-6.67707 


150.385 


-623.767 


858. / / 7 


7 


4.90673278031246 


11.9835 


-129.523 


757.114 


-1380.49 














oo 


4 











Table 3.7: 1/N expansion coefficients for 14 when F = 8. 
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N 


76= a + b/N + c/N'^ + d/N^ 






a 


b 


c 


d 


3 


103.537190082645 










4 


57.4074000134889 










5 


30.083006000112 










6 


8.52435062302512 


-202.192 


2030.17 


-5851.51 


7537.7 


7 


-1.73452473741117 


248.055 


-4723.54 


27466.8 


-4649 














oo 


-24 











Table 3.8: 1/N expansion coefficients for Iq when F = 8. 



N 


Jg = a + b/N + c/N^ + d/A^=^ 






a 


b 


c 


d 


3 


1073.4573622182 










4 


-228.646302651363 










5 


-1293.26989559543 










6 


-1868.87943727765 


1402.47 


-68534 


384230 


-544767 


7 


-1927.54995720038 


13256.9 


246350 


1.26146 X 10^ 


-1.9673 X 10^ 



Table 3.9: 1/A^ expansion coefficients for /g when F = 8. 
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4 Diagrammatic expansions 



The proof of the sixth moment sum rule [13] is based on an analysis of the Mayer diagrammatic 
expansion of the direct correlation function, in the fiat space in the bulk. Therefore, it is 
interesting to study this expansion on the sphere. 

The direct correlation function c is defined by the Ornstein-Zernike equation 

h{ri, fs) = c(fi, fs) + p j h{fi, r3)c(f3, f2)dS3, (4.1) 

where fi points the direction of a point in the sphere from its center. Since the system is 
homogeneous, both h{fi,f2) and c(ri,r2) depend only the angle 9i2 = (^i,'r2). It is useful to 
introduce the expansion of h and c in Legendre polynomials 

oo 

h{9) = J2hiPiicose) (4.2) 

1=0 

hi = '^^ j h{x)Pi{x)dx (4.3) 
and similarly for c. The Ornstein-Zernike equation reads 

hi = ci + -^p^hici . (4.4) 

implying 

hi = \— . (4.5) 

Notice also that the moments l2n of the total correlation function are directly related to the 
coefficients hi, since P;(cos^) can be expressed as a polynomial in sin^(^/2) using the formula 
cos^ = 1 -2sin2(^/2). 

As explained in detail in [13], the direct correlation function has a renormalized diagram- 
matic expansion which reads 

oo 

c(ri,r2) = -Tv{ri,r2) + c^°\ri,r2) + ^ c'^'*''(n, ^2) (4.6) 

s=l 

where v{fi,r2) is the Coulomb potential between to unit charges. c*^°^ is the renormalized 
"watermelon" Meeron graph 

c^''\fi,f2) = lK\f,,f2) (4.7) 
where K is defined by the equation 

K{ri,f2) = -Tv{fi,f2) - p j Tv{fi,r3)K{f3,fi)dS3. (4.8) 

The rest of the expansion of c, Xl^i c^^\ri,r2) correspond to the remaining diagrams; see [T3] 
for details. 
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The coefficients Ki of the expansion of K in Legendre polynomials can be obtained from (14 .Sp 



where the coefficients of the Coulomb potential are obtained by solving Poisson equation in the 
sphere 

2/ + 1 , , , 

Notice that the coefficient for / = is not defined properly, rather as shown in [20j it is a 
constant that can be chosen arbitrary. Then 

In [13] it is proved that the zeroth and second moment of each c*-^-* vanishes, therefore they 
do not contribute to small-/c expansion of the Fourier transform of c up to terms of order /c^. 
Then the structure of the Ornstein-Zernike equation in Fourier space guaranties that these 
terms do not contribute to the expansion of the total correlation function h up to k^. This fixes 
the value of the moments up to the sixth. Unfortunately, on the sphere this argument breaks 
down. If one can show that the zeroth and second moment of c^^^ also vanish on the sphere, 
this fixes the values of cq, and ci, and from eq. (14.41) one can deduce the value of ho, hi and 
therefore /q and I2 but unfortunately one cannot deduce anything exact regarding /i2, h^, and 

4.1 Second moment 

Suppose, as in the plane, that c^'^^ has vanishing zeroth and second moment, i.e. Cq^ = 0, and 
c^^ = 0. Then to compute q, eq. (14. 6 p with the sum ignored implies 

ci = -Tvi + cf = -Vvi + ^mi . (4.12) 

The coefficients c[°^ = mi/2 are obtained from those of the expansion of K"^ , 

00 

K'^{e) = ^miPi{cose) (4.13) 
1=0 

which can be obtained if the expansion of a product of two Legendre functions is known. If 



i"=0 

with 



Pi{x)Pv{x) = Y,P%Pi"{^) (4.14) 
P\'' = ^^j Pi"{x)Pu{x)Pi{x)dx (4.15) 
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then 



1=0 l'=0 



The p['i', can be computed [llj using the recursion equation 

{21 + l)xPi{x) = IPi-iix) + (/ + l)Pi+i{x) . 
For instance, for /" = 1, using x = Pi{x) in 04.171] and replacing it in fl4.15l) gives [TT] 



2(/ + r 



2 [(2/ + l)(2/ + 3) 



21 



{21 - 1)(2/ + 1) 



Therefore, 



mi 



1=0 



(2/ + l)(2/ + 3) 



/ + 1 



^ {2{l + 1)/ + NT){2{1 + 1)(/ + 2) + A^r) 



2 

3r 

2N' 



1=0 



1 1 
_2l{l + l) + NT ~ 2{l + l){l + 2) + NT_ 



Then, replacing ( KW\ . into f l412l) and (iS 



ar 
















r 


1 + 


4 


~ 4 



But 



hi = ^ j h{e){l-2sm\e/2))d{cosi 



(4.16) 



(4.17) 



(4.18) 



(4.19) 



(4.20) 



(4.21) 



Since Iq = —1, we verify that fl4.20p is compatible with the sum rule from section ^ I2 = 
NT/{T - NT -A). 

Since the second moment sum rule was proven on firm grounds in section [31 going backwards, 
this proves indirectly that '^i'''' = 0. 
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4.2 Higher moments 

Although there is no reason to suppose that the contributions from Yl^i 4*'* for / = 2 and / = 3 
will vanish for a finite sphere, it is interesting to compute the contribution from —Tvi + c^^^ to 
the 4th and 6th moments. Let us define an approximate direct correlation function c'^pp''^^ by 



(4.22) 



/=0 



with 



^approx 



-Tvi + -mi 



(4.23) 



and the corresponding total correlation function /^^pp™^ obtained from Ornstein-Zernike equa- 
tion 

-(2/+l)r mi 



^approx 



1 + 



Nmi 

2/(/ + 1) ~ 2(2/ + 1) 



4.3 Fourth moment 

The fourth moment is related to the coefficient h2 by 



and so 



5pb 
N 
5 
N 



h{d) 



l-^(2sin(0/2))2 + ^(2sin(e/2))^ 



' AT + (ATF 



(A^r)^ 

24 



+ 1 + 



12 



r - A^r - 4 



The quantity h2^^^°^ is given by 



dS 



(4.25) 



(4.26) 



^approx 



Sr _|_ m2 

12 ' 2 

~'~ 12 10 



with 



and HI 



m2 = J2J2^^^^'P2 



1=0 l'=0 



6(/ + l)(/ + 2) 



(2/+ l)(2/ + 3)(2/ + 5) 



Sl',l+2 + 



6(r + i)(r + 2) 



(2/' + l)(2/' + 3)(2/' + 5) 



SlI'+2 — 



(4.27) 



(4.28) 



2/ + 1 



. (4.29) 
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Then 



with 



1=0 



7712 



/=0 



(2/+ l)(2/ + 3)(2/ + 5) 



+E 

1=0 



5Kfl{l + 1) 



(2/-l)(2/ + l)(2/ + 3) 



(4.30) 



15(/ + !)(/ + 2) 



+ 



5(2/ + l)/(/ + l) 



(2/(/ + 1) + x)(2(/ + 2)(/ + 3) + x)(2/ + 3) (2/(/ + 1) + s)2(2/ - 1)(2/ + 3)_ 

(4.31) 



Replacing this into f l4.27p and fl4.26p gives an approximation for the fourth moment /, 



approx 
4 



4.4 Sixth moment 

Following similar steps, one can obtain an approximation for Jg. The sixth moment is related 
to hs by 



f 24/2 120/4 160/6 



Therefore, using f l4.26p and the known values of Iq and 12, 



(NTf 



Nh Nh2 1 



9 



1120 160 40 40(r-A^r-4) 



h2 can be approximated by /ig ^ from last section, and can be approximated by 



^approx 



24 2 

1 _i_ EL _ A^"i3 

24 14 



with 



and 



1=0 l'=0 



(5x3 -3x)P/(x)Pr(x) 



(4.32) 



(4.33) 



(4.34) 



(4.35) 



(4.36) 



P3' can be computed by applying the recursion equation (14.1 7p several times to obtain 

1 



21 + 1 



{21 - 1){21 - 3)' 



(/ + l)(/ + 2)(/ + 3) 
(2/ + 3)(2/ + 5) 



^^+3(2;) + 



(2/-3)(2/ + 3) 
3(/ + l)((/ + 1)^-2) 
(2/- l)(2/ + 5) 



(4.37) 
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Then replacing into (14.361) 



7 



2(2/ + 1) 



5/(/ -!)(/- 2) 
(2/- l)(2/-3)(2/-5) 



3/(/2 



5(/ + l)(Z + 2)(/ + 3) ^ 3(/ + !)((/ + 1)2 



2/ - 3)(2/ + 3)(2/ 
1) 



(2/ + 3)(2/ + 5)(2/ + 7) 



(2/- l)(2/ + 5)(2/ + 3) 



6r 



Replacing this result into (I4.35P one obtains 

(/ + !)(/ + 2) 



1=0 



31 



{21 + 5)(2/(/ + 1) + iVr) [(2(/ + 1)(/ + 2) + iVr)(2/ 
5(/ + 3) 



(2(/ + 3)(/ + 4) + Arr)(2/ + 3) 

7(2 + 3A^r)r 



2A^(12 + 3A^r) 



Then 



, approx ^ -7(-4 + iv(4-6r) + iv2r)r 

^ A^(96 + 4(7A^ - l)r + (6 - Ar)Arr2) 



(4.38) 



(4.39) 



(4.40) 



^approx /i^PP"^"^ can be replaced into f l4.33p to obtain an approximation for the sixth moment 



rapprox 
-'6 



4.5 Comparison with numerical results 

The expressions for the fourth and sixth moments can be compared to their exact values which 
have been computed numerically in section 13.21 Tables 14. 14. 2[ 14. 3[ 14.41 show the exact value 
of the moments and the approximate value obtained from the diagrammatic technique. For 
the fourth moment j^pp'""^ gives an acceptable approximation, on the other hand j^pp''"^ has 
accumulated too many errors from h^'^^'^^ and h^^^^^ and gives a poor agreement with the 
exact value. However, as iV — )■ cxd both /|'PP''°^ and /^pp''"^ approach the known thermodynamic 
limit, as expected. 
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h 


7-api)rox 
i 


I i error 


h 


7-api)iox 
-'() 


/(i error 
















2 


-1.06666666666667 


-0.686993 


35.6% 





0.241677 




3 


-0.73469387755102 


-0.534741 


27.2% 


3.30612244897959 


1.54027 


53.4% 


4 


-0.552915766738661 


-0.42966 


22.3% 


5.1605471562275 


2.52407 


51.1% 


5 


-0.437781621713968 


-0.356903 


18.5% 


6.11292630320115 


3.21788 


47.4% 


6 


-0.361584090880502 


-0.304426 


15.8% 


6.56602660174796 


3.71182 


43.5% 


7 


-0.307439760694233 


-0.26506 


13.8% 


6.78048920359037 


4.07309 


39.9% 


8 


-0.267158562552772 


-0.234542 


12.2% 


6.87498575554362 


4.34502 


36.8% 


9 


-0.236094785664912 


-0.210236 


11.0% 


6.9075123159669 


4.5552 


34.1% 


10 


-0.211435346122364 


-0.190444 


9.9% 


6.90772192712509 


4.72147 


31.6% 


11 


-0.191399568920847 


-0.174026 


9.1% 


6.89103601184945 


4.85568 


29.5% 


12 


-0.17480728582518 


-0.160195 


8.4% 


6.86574608931141 


4.96591 


27.7% 


13 


-0.160846001528985 


-0.148387 


7.7% 


6.83643498600859 


5.05783 


26.0% 


14 


-0.148938895568825 


-0.138191 


7.2% 


6.80567000027031 


5.1355 


24.5% 
















oo 
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Table 4.1: Exact value of the moments and their approximation from the diagrammatic tech- 
nique for r = 4. 



N 


TT. 

h 


fapprox 
-'4 


li error 


h 


fapprox 
-'6 


— ^ 

Iq error 
















3 


1.5 


2.12597868844099 


41.7% 


32.4 


7.00428372252487 


78.4% 


4 


1.77112299465241 


2.11737368908595 


19.5% 


29.2857260386672 


5.04095606861884 


82.8% 


5 


1.92727455514225 


2.10358894768327 


9.1% 


24.3553395496129 


3.85885436737017 


84.2% 


6 


1.97167595506035 


2.09127008496943 


6.1% 


19.3570221116088 


3.09958641182161 


84.0% 


7 


1.99917180579183 


2.08108076102509 


4.1% 


15.8157978386411 


2.57924397229046 


83.7% 


8 


2.01319334170647 


2.07273170306951 


3.0% 


13.2137920826172 


2.20346441160235 


83.3% 


9 


2.02020128573558 


2.06584181515105 


2.3% 


11.2455089383337 


1.92065635218756 


82.9% 


10 


2.02414460391682 


2.06009120107809 


1.8% 


9.75173060605997 


1.70073355331074 


82.6% 


11 


2.02625080279343 


2.0552339688013 


1.4% 


8.58562698480205 


1.52513914797188 


82.2% 


12 


2.02719976576841 


2.05108476391051 


1.2% 


7.65247243254821 


1.38187395446424 


81.9% 
















oo 


2 














Table 4.2: Exact value of the moments and their approximation from the diagrammatic tech- 
nique for r = 6. 
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N 


h 


7-approx 
-'4 


I4 error 


h 


T-approx 
-'6 


Iq error 
















3 


5.9504132231405 


7.336835008 


23.3% 


103.537190082645 


-127.1369698 


222.8% 


1 


5.3521(3711227873 


6.1509183 


11.9'X 


57.1071000131889 


-71.80581066 


225.1% 


5 


5.31949584221798 


5.585392260 


5.0% 


30.083006000112 


-54.66045440 


281.7% 


6 


5.03608123946663 


5.254867561 


4.3% 


8.52435062302512 


-46.46113141 


645.0% 


7 


4.90673278031246 


5.038212233 


2.7% 


-1.73452473741117 


-41.68739389 


2303.4% 


8 




4.885278337 






-38.57330645 




9 




4.771577915 






-36.38521278 




10 




4.683739174 






-34.76512324 




11 




4.613844211 






-33.51802809 




12 




4.556906725 






-32.52880274 




13 




4.509630618 






-31.72517539 




14 




4.469749993 






-31.05951906 




15 




4.435655939 






-30.49918829 




16 




4.406174568 






-30.02106637 




17 




4.380429327 






-29.60832983 




18 




4.357752484 






-29.24844597 




19 




4.337626449 






-28.93188851 




20 




4.319643825 






-28.65128805 


















oo 


4 
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Table 4.3: Exact value of the moments and their approximation from the diagrammatic tech- 
nique for r = 8. 
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N 


h 


rapprox 
-'4 


I4 error 




7-approx 
-'6 


Iq error 




2 






3 




















2 


-0.666666666666667 


-0.6317574181 


5.236% 


-0.8 


-0.752415111 


5.95% 


3 


-0.9 


-0.8700339821 


3.330% 


-1.35 


-1.337880192 


0.90% 


4 


-1.06666666666667 


-1.042131679 


2.300% 


-1.82857142857143 


-1.864218645 


1.95% 


5 


-1.19047619047619 


-1.170437198 


1.683% 


-2.23214285714286 


-2.312761744 


3.61% 


6 


-1.28571428571429 


-1.26919496 


1.285% 


-2.57142857142857 


-2.689601017 


4.60% 


7 


-1.36111111111111 


-1.347327155 


1.013% 


-2.85833333333333 


-3.00618709 


5.17% 


8 


-1.42222222222222 


-1.410579495 


0.819% 


-3.1030303030303 


-3.27362061 


5.50% 


9 


-1.17272727272727 


-1.162780191 


0.675% 


-3.31363636363636 


-3.501250962 


5.66%) 


10 


-1.51515151515152 


-1.50656529 


0.567% 


-3.4965034965035 


-3.69658924 


5.72% 


11 


-1.55128205128205 


-1.543801073 


0.482% 


-3.65659340659341 


-3.865574853 


5.72% 


12 


-1.58241758241758 


-1.57584505 


0.415% 


-3.7978021978022 


-4.012891553 


5.66% 


13 


-1.60952380952381 


-1.603706117 


0.361% 


-3.92321428571429 


-4.142245077 


5.58% 


14 


-1.63333333333333 


-1.628149088 


0.317% 


-4.03529411764706 


-4.256585403 


5.48% 


15 


-1.65441176470588 


-1.649763932 


0.281% 


-4.13602941176471 


-4.358278405 


5.37% 


16 


-1.67320261437909 


-1.669012791 


0.250% 


-4.22703818369453 


-4.449236580 


5.26% 


17 


-1.69005847953216 


-1.686262682 


0.225% 


-4.30964912280702 


-4.531018304 


5.14% 


18 


-1.70526315789474 


-1.701808689 


0.203% 


-4.38496240601504 


-4.604903343 


5.02% 


19 


-1.71904761904762 


-1.715890707 


0.184% 


-4.4538961038961 


-4.671950573 


4.90% 


20 


-1.73160173160173 


-1.72870573 


0.167% 


-4.51722190852626 


-4.733042400 


4.78% 


21 


-1.74308300395257 


-1.740417020 


0.153% 


-4.57559288537549 


-4.788919165 


4.66% 


22 


-1.7536231884058 


-1.751160988 


0.140% 


-4.6295652173913 


-4.840206004 


4.55% 


23 


-1.76333333333333 


-1.761052508 


0.129% 


-4.67961538461538 


-4.88743397 


4.44% 


24 


-1.77230769230769 


-1.770188989 


0.120% 


-4.72615384615385 


-4.931056780 


4.34% 


25 


-1.78062678062678 


-1.778653568 


0.111% 


-4.76953601953602 


-4.971464120 


4.23% 


26 


-1.78835978835979 


-1.786517625 


0.103% 


-1.81007115189871 


-5.008992111 


1.11% 


27 


-1.79556650246305 


-1.793842793 


0.096% 


-4.84802955665025 


-5.043933457 


4.04% 


28 


-1.80229885057471 


-1.800682557 


0.090% 


-4.88364849833148 


-5.076541501 


3.95% 


29 


-1.80860215053763 


-1.807083561 


0.084% 


-4.91713709677419 


-5.107039000 


3.86% 


30 


-1.81451612903226 


-1.813086665 


0.079% 


-4.94868035190616 


-5.135621381 


3.78% 


31 


-1.82007575757576 


-1.818727815 


0.074% 


-4.97844251336898 


-5.16246097 


3.70% 


32 


-1.825311942959 


-1.824038756 


0.070% 


-5.00656990068755 


-5.187710286 


3.62% 
















oo 


-2 






-6 







Table 4.4: Exact value of the moments and their approximation from the diagrammatic tech- 
nique for r = 2. 
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5 Fluctuation of linear statistics 



It was remarked in the Introduction that the sum rules for the moments of p^-^ (f, 0) are related 

to the average value of the linear statistic Ylh=i l^d^"- Here we point out that in the case of the 
soft disk geometry and n = 1 the exact distribution function for this linear statistic is easily 
evaluated. Two features of the exact distribution are singled out for further study. One relates 
to the 0(1) correction term to the mean. We are able to deduce the value of this correction 
to the mean for a general smooth radial linear statistic. The other is that for large the 
exact distribution tends to a Gaussian, with mean and variance as can be anticipated from a 
Coulomb gas viewpoint [8]. We will show that a mechanism for such Gaussian fluctuations in 
the case of general rotationally invariant linear statistics can be identified from the viewpoint 
of the expansions (12.31) and (12.81) . 

First, with (■)piane denoting an average with respect to the Boltzmann factor (II. Sp . we 
observe that a simple scaling of the coordinates fj gives 

/ plane V LTTph/ 

Coulomb gas theory ^ tells us that setting R = 1, so pb = N/n, then taking iV — )■ oo, the 
distribution function of a radial linear statistic A = Y^^=i9i''^i) should become the Gaussian 

with 



plane N—^oo 



Mn = 2N [ rg{r)dr + 0(1) 
Jo 



'0 

1 



^2 = - / r{g'{r)fdr + o{l). (5.3) 





5.1 Universal correction to the mean A4n 

According to (15. ip . with = N/tt, 

This is of the form (15. 2p with 

a' = i (5.6) 

and (15. 5p . (15.60 are in turn consistent with (15. 3p in the case g{r) = r^. The expression (15.50 is 
an identity for all N. But with g{r) = r^. 



Mn= / r>(i)(r)rfV. (5.7) 

JR2 
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Writing = PbXo<r<R, where xa = 1 if ^ is true, and Xa = otherwise, this can be written 

Mn = ^R^Pb + [ r=^(P(i)(r1 - Pb) (fr. (5.8) 

Recalhng that in (15. Sp . pb = N/n and R = 1, comparison of (15.80 and (15. 5p gives the sum rule 

J^fiPiim-Pb)d'f= - (5.9) 

We remark that an alternative derivation of (15. 9p . valid for F even, is to use (12. 4p and (12. 9p . 
together with the recurrence relation for the gamma function. We remark too that (15. 9p is the 
soft disk version of the sum rule for the 2dOCP in a disk with hard wall boundary conditions 

P(i)(0)-(l-^)p, = -^^ /" r'(^p^,^{R-r)-p,)dr, (5.10) 



where P{i){r) is measured inward from the boundary. This latter viewpoint allows the N ^ oo 
limit of (15. 9p to be formulated. Thus on the LHS change to polar coordinates, and further 
change variables r h-> i? — r. Defining P(7)(^) ~ P(i)(-^ ~ ''^); corresponding to the density as 
measured inward from r = i? in the soft wall disk OCP, and noting that charge neutrality 
requires 

[ (R-r)ip^^^(r)-p~,)dr = (5.11) 

we see that 

oo pR 

r3(p(i)(r) - p,) dr= (R - rfip^^^ir) - p,) dr 

J — oo 

r-R 



I (-2i?V + 3i?r2 - r3)(p^7)(r) - p^) dr 

J — oo 

/oo 
rip'(J)ir) - Pb) dr + OiR) (5.12) 
-oo 



(compare with working below (3.17) of [I9]). Comparison with (15. 9p we thus have that in the 
thermodynamic limit 

-2nT r r(p^7)(r) - p,) rfr = (l - (5.13) 

J — oo ^ 

The sum rule (I5.13P should be compared against the contact theorem [7] for the hard wall 
plasma, 

pS(0) - 27rr r(pfj^(r) - p,) rfr = (l-^), (5.14) 

where P('^(r) denotes the density as measured from the boundary. 

The above working suggests that the first two terms in the asymptotic expansion of A^tv in 
the case of g{r) = r™, for m G Z+, admit a form analogous to that exhibited in (15. 5p . First, 
with this choice of g{r), analogous to (15. 8 p we have 

o poo 

Mn = -—-R'^-'^Pb + 2n / r"^+^(p(i)(r) - pb) dr . 
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Manipulation of the integral as in f l5.12p now shows 

o poo 



Substituting of (I5.13P we then obtain 



= J^R-'^'P^ + ^(l - ^R"" + 0{R-'pl"\ (5.15) 
In particular, with R = 1 and pb = N/tt, we have 

In this setting write P(i)(r) = {N/iT)xQ<r<i + f^i^) + 0{N^^^'^). Then it follows by comparing 
the definition of A^at in the case g{r) = r"^ to f l5.16p that 



27r j r"'+^K{r)dr = -(l-- 



And it follows from this that K{r) is the distribution supported at r = 1 given by 

Thus for g{r) smooth in the neighbourhood of r = 1, and with R = 1 and pt = N/tt, we have 
that for large 



Mn = 2nN rg{r) dr+^(l- ^) (^'(l) + 0(Ar-i/2). 



(5.18) 



5.2 Derivation using ( 12.91) 

We will now consider the structure of the Gaussian fluctuation formula fl5.2p itself. It can 
readily be deduced at F = 2 as a consequence of the expansion formula (12. 9p for lN,r[g] [H]- 
Thus then p = 1 and cl^\l) = 1 for /i = 6n, and clf^\l) = otherwise. It is interesting to 
probe how (12. 9p can lead to (15. 2p for general F = 4p + 2. We then have 



plane /7v,rilJ 



where, with specified by (12. 5p 

2 ^ 

/^_p[e^MO] = iVIvr^J^ (c(f)(2p+ 1)) J] ^'^Je-^^^'/'e^'^^^^')], (5.20) 

(U 1 = 1 

and according to (12. 6p 

^(A^-l)>/ii>/i2>--->/iiv>0 (5.21) 
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with \fi\ = rN{N — l)/4. Let's introduce the scaled variables fti := 2^i/NT so that to leading 
order in N 

1 > /il > /i2 > • • ■ > /iAf > 0. 

In terms of these scaled variables, and upon the change of variables = s, we see 



oo 







For large the maximum of the exponent occurs aX s = jli. Expanding the integrand about 
this point and completing the square shows that for large 

G,[e-^^^V2,.M.^)] ^ G,Je-^^^^/^]exp (^ka{^^ - ^{a{^i)f). (5-22) 



4rA^ 

To proceed further we hypothesize that the values of {pn} for which the cl!^\2p + 1) in (12.91) 
are nonzero are uniformly distributed (or at least this situation dominates). Then we have 

N 

(=1 

~ n G/., [e-^^"'/2] exp (ikN / a(v^) du - -- (a'(v^))2 du) . (5.23) 
Substituting in (I5.20p then substituting the result in (I5.19P reclaims (15. 2p . 
5.3 Test near T = 2 

In this subsection we verify (I5.16P for F = 2 and for F close to 2, where explicit analytical 
calculations of A^at can be done. When F = 2, the density profile is [9l §15] 

,,.4, 

where upb = N. Then, with m = 2n, 

f 2n , ^ 1 ^ {k + n)\ N{N + n)\ 
Mn = / P(i)(r)d^r = > = ) -—- . (5.25) 



Now, when N ^ oo, 

, N n n(n — l)(3n + 2) ^ , , , s 

MN\r=2 = + 2 + ^ ^ + ^^-26) 

which satisfies (I5.16P when F = 2. It is interesting to notice that there are no 0(A^~^/^) 
corrections to A^at at F = 2. This is, however, a special feature of the F — Z C3iSG, clS for general 
F there will be non zero 0(A^~^/^) corrections (see sec. 15.41) . A similar situation happens for 
the finite-size corrections to the free energy of the 2dOCP in the soft disk at F = 2 [19j . 
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For r close to 2, one can perform an expansion in powers of F — 2 of the density profile, to 
compute AiN, similar to a computation done in [12j. To the linear order in F — 2, the density 
is given by P(i)(r; F) = P(i)(r; 2) — ^^{p{r)(3H)'^ + 0((F — 2)^), where the truncated average is 
taken with the Boltzmann factor f ll.Sp at F = 2, i7 is the potential energy of the 2dOCP and 
p(r) is the microscopic density. Explicitly, 



R2xR2 



[p(3)(fi,f2,f) - P(i)(f)p(2)(fi,f2)] f (fi,f2) (ffid^f2 + 7rpferV(i) (r) 



+ 2/ p^2)in,r)v{fi,r)(ffi + npb [p(2)(fi, f) - p(i)(fi)p(i)(f)] 



(5.27) 

where the density and correlations in the RHS are evaluated at F = 2 and v{f, r ') is the Coulomb 
pair potential (11.11) . At F = 2, the correlation functions have the simple structure [12] [9], §15.3] 



P(i){ri,r2, ...,re)= pldet{K{zi, Zj))i<ij<e , 



(5.28) 



with 



N-l 



k=0 



(^) 



(5.29) 



and z = y/Wpbre^'^, where (r, 0) are the polar coordinates of r. To perform the integrals in (I5.27p . 
it is convenient to expand the Coulomb potential (11.11) in a Fourier series in the polar angle 



°° 1 /r A"" 
;(fi, r2) = — — cos[m(0i - ^2)] - log ■ 

m=l \ 



(5.30) 



with r> = max(ri,r2) and r< = min(ri,r2). Replacing (15.281) and (I5.30p in (I5.27p . after much 
algebra, one obtains 



P(i)(^;r) = P(i)(f;2) 



N-l N-l 



\2k2 



N-l N-l 



1=0 fc2=0 J- V ^ / fc^=ofc2=fcl+l ^ V ^ 



I^|2fc2 l^p'^r 



+ 



fe2#fci 



^-1 |^|2fc 



fcl=0 fc2=fci+l 



17(^2 + 1, |zn + |z|2'=2F(A;i + 1,1^1 



(^2 - A;i)A;i!fc2! 



TV-l 



^ 2fco! ^ 



2fc2 ^1 7(A;i, logd^n + Jj J e-*t^i logt rft 



N-l 



|2fei 



fc2=0 
Af-1 



fci=0 



fci=0 



-2 



2A,i 



2(A;i!) 



fei=0 

where we have defined 



7(fci + l, |2niog(|zn+ / e-H^'\ogtdt 



I{ki,k2)= [[ e-'^-'H'lH'^^ dtidt2 

J Jo<t2<ti 
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+ 0((F-2)^ 



(5.31) 



(5.32) 



and 



noo 
e-''-H\Hl^ log(max(ti, ts)) dtidt2 . (5.33) 

Some useful properties, as well as the asymptotic expansion when /ci — )■ oo and k2 oo oi 
I{ki, ^2) and J{ki, ^2), are discussed in the appendix. 

Consequently, with Tcp^ = N, the 2n- moment of the density profile is 



M 



I fci=0 ^' 

^ ^1 ^1 (^x(A;^,fc2)_x(A;i + n,A;2) + ^X(fci,fc2)-X(fci,A;2 + n) 

ki=0 fc2=fcl + l 

N-l N-1 



h\k2\{k2 - k^ 



+ EE 



2ki\k2\ 



{ki+n)\ 
kil 



J{kuk2) - J{ki + n, k2) 



+ 0((r-2)2). (5.34) 



In the last term of ( 15.34^ it is convenient to use the recurrence relation (]A.6l) and write 



{ki + n)\ 



rt-l 



J{ki,k2)- J{ki + n,k2) 



V (fci+n)! 
{ki + n)\ 

{kl + l + ^)\ 



n-l 

E 

£=0 



x{k, + iM)- ^^'^P' AkiM) 



k,\ 



(5.35) 



Consequently A^at is given by A^jv = Mn\y=2 ^ (^ — 2) (mi + m2 + ms) + 0((r — 2)^) with 



AT-l N-l N-l n-l s. 



A;i=0 



ki\k2\ ^ A;i!(fci + ! + £)' 



(5.36) 



m2 



■ 2 ^ ^ fcilfcs! ^ (A;i + 1 + 

A;i=0 fc2=0 £=0 ^ 



X(fci+^,fc2)- ^V/^' ^(^l'^2) 



(5.37) 



fe2#fei 



and 



^1 ^1 ^X(fci, A;2) - X(fci + n, ^2) + ^X(fci, ^2) - X(A;i, ^2 + ri) 

fci=0 A;2=A;i + l 



kx\k2\{k2-ki 



(5.38) 

In f l5.36p it is convenient to split the sum on k\ and /c2 into two sums for k\ < k2 and ki > k2 
and in the latter use (1A.3I1 to find rhi = + m\ with 



fci=0 



n-l 

*.E 

i=0 



jki+n)] 



(5.39) 
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and 



N-^^ ^1 j^^^^^^^nz^/ (^k^ + n)\ {k2 + n)\ 



In (15 .39 p . the expression inside the sum is a polynomial of ki of order n — 1, therefore to 
obtain the order 0{N^) when — )■ oo of it is sufficient to consider the leading order of the 
polynomial and replace the sum over ki by an integral, 



AT-n j-N -ri-l 

- / dk^Y,[{^ + l)k^^-' + o{k'i-^ 



m 



2 

n + l 



4 

In the appendix it is shown that, for large ki and ^2, 

X(fci,/c2) 1 ^ / k2-ki 



+ 0(1/Ar) . (5.41) 



ki\k2\ 2 \ ^2(A;i + A;2) 



(5.42) 



To obtain the leading order contribution in of fh\ one should replace the sums over ki and 
^2 in (15.401) by integrals, and realise from (I5.42p that the leading contribution will be obtained 
for ki and ^2 laying in a strip following the line ki = k2, with ki < k2, and of width of order 
'IV. With this in mind, it is useful to do a change of variables in the integrals, from ki and /c2 



to k^ = ki + ^2, which goes along the strip, and u = {k2 — ki)^j2k^, which goes perpendicular 
to the strip. The integral over u will give 0(1) contributions, due to the fast decay of the erfc 
function, while the one over fc+ will give a contribution of order 0(iV"). To obtain this, it is 
sufficient to keep the leading order in fc+ in the integrand. In (I5.40p . the expression in the sum 
over £ is a polynomial of ki and k2 that vanishes when ki = k2, therefore it can be written as 



with Q{ki, /C2) a polynomial in ki and k2 of order n — 2. Actually, to the leading order. 



n-2 



^2) = n k^ ^ '7^2 + polynomial in ki and k2 of order n — 3 



r=0 

n-2 



^ 2 

Then, replacing in (I5.40p . 



n(n - 1) ( ^ ) + 0{kl-^) . (5.44) 



N- 



2 ^ ^ ki\k2\ 

ki=0k2=ki+l 
r-oo r2N 



4 

n — 1 



du dk+ueiic{u) [n{n - 1)2-^'^-^'^ k^'^ + Oikl'^)] + 0(1/ ViV) 
Jo Jo 

+ 0(l/yiV), (5.45) 



30 



where the integral ueifc^u) du = 1/4 was used. Putting together the last results, we find 

mi = ^ + 0(l/yiV) . (5.46) 

Regarding the other contributions to AiN, it can be noticed using the recurrence rela- 
tion (lA.Sp that m2 + rfi^ = when n = 0, 1, 2, and 3. For n > 3, the leading order of these 
terms can be obtained following similar steps as the ones done for rhi. The sums over ki and 
k2 are replaced by integrals over the variables /c+ and u defined above, and only the leading 
order in /c+ is kept. From flA.131) . one can obtain that, for fc+ large. 



k2\{k, + e)\ 



Substituting into (15.371) shows 



rh2 = -iV-"2 / du 

J-oo Jo 

n — 1 



2N 



dk-L. \ l 2/u_i 



/e"" (k 



n-l 



+ 0(1/VA^) 



o(iA+) 



n—l 



-±] +0(l/v^) 



(5.47) 



(5.48) 



For m3, use 



(^X(fci, k2) - I{k, + n, k2) _ klne-' 



ki\k2\ 



2"y2^ 



and 



I{ki,k2) -I{ki,k2 + n) fc"ne-"' 



ki\k2 



2'^y2^ 



-i + ^^^ + o(iA, 



i + ^^ + o(i/M 



to obtain 



A^-" / du 
n—l 







2N 



dk+ a/2A;+ f k^\^ 2e ^'^n{n — l)u 
u 



Oil/VN) 



4 



+ Oil/VN). 



(5.49) 



(5.50) 



(5.51) 



Therefore m2 + rfis = 0{l/yN) do not contribute to order 0(1) in A^a?. Summing up all 
results, 



M 



N 



N 



^ ^ ^ + - - (F - 2)- + 0(l/v^) + 0((F - 2f) 



(5.52) 



This result is in agreement with the general formula fl5.16p when it is expanded around F = 2 
to the linear order in F — 2 with m = 2n. 
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5.4 Numerical results for F = 4, 6, and 8 



For r even, by application of f l2.4p or fl2.9p . the 2?T,-moment of the density can be expressed as 



M 



N 



{NT/2y 



^soft 



AN) 



(r/2)) 



2 N 

e=i k 



N 



(/ife + n)\ 



with 



^soft 



V (cr(r/2))^ fr ■ 



(5.53) 



(5.54) 



which is the partition function of the 2dOCP in the soft disk, up to a muhiphcative constant. 
Tables [5m - l5. 91 show the numerical evaluation of at as well fit to Mn = aN + b+cN-^/^ + 
d/N, for r = 4, 6, and 8, and n = 2, 3, and 4. These numerical results, and the proposed fit, 
show that fl5.16p is indeed satisfied. The convergence is very good for all moments for F = 4 
and 6. For F = 8, unfortunately we were not able to perform the calculations beyond N = 7 
particles, the results presented have not yet converged to the N ^ oo expected value. Notice 
that, contrary to the case F = 2, there are non-zero 0(A^~^/^) finite-size corrections to Aij\i- 



N 



M 



N 



aN + 



c/^/N + 



d/N 







a 


h 


c 


d 


2 


0.8125 










3 


1.126262 










4 


1.44563743218807 










5 


1.76891109591098 


0.336437 


-0.0742714 


0.458953 


-0.221263 


6 


2.09454890418255 


0.333419 


-0.00242266 


0.269645 


-0.0817588 


7 


2.42171814295119 


0.332974 


0.0108746 


0.230301 


-0.0491322 


8 


2.74996856529295 


0.33338 


-0.003688474 


0.277614 


-0.0922608 


9 


3.07901955876735 


0.333458 


-0.00698637 


0.289203 


-0.103695 


10 


3.40868118671838 


0.333396 


-0.00397741 


0.277889 


-0.0917455 


11 


3.73882025776555 


0.333353 


-0.00169446 


0.268779 


-0.0815306 


12 


4.06934089384864 


0.333341 


-0.000983289 


0.265786 


-0.0779911 


13 


4.4001722794716 


0.333341 


-0.000982974 


0.265784 


-0.0779894 


14 


4.73126081887937 


0.333342 


-0.00105077 


0.266097 


-0.0783948 














oo 




1/3 










Table 5.1: Fourth moment (n = 2) of the density when F = 4. 
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N 


Mn = aTV + b + c/y/N + d/iV 






Qj 


u 


C 




2 












3 


l.UOOOOOOOOOOOOO 










4 












5 


1 59318930981443 




04991 87 


\J .OOOiyO 


-0 07451 99 


6 










n 070(^77'^ 


7 


1 Q^Q1R9'^^nn'^44Q 

J- . z/OiJ ±\J^OO\J\JO^^z/ 




n 0?^*^! 7^R 






8 


9 99RRnni A791 04^ 










9 


9 4(S^Q^1 Q79^7 




n 09199"^ 




n 91 Q4nR 


10 


9 7n(S7(^'^nRzL9^79zL 


D 9zLQ7zl^ 






n 991 "^4^ 


11 


2.94871984780432 


0.249774 


0.0193389 


0.671799 


-0.227593 


12 


3.19159601492947 


0.249824 


0.0163325 


0.684451 


-0.242556 


13 


3.43522469549061 


0.249871 


0.0132299 


0.69815 


-0.259556 


14 


3.67947926593368 


0.249905 


0.0108374 


0.709186 


-0.273865 














oo 




0.25 









Table 5.2: Sixth moment {n — 3) of the density when F = 4. 




Mn = aN + b + c/VN + d/N 






a 


b 


c 


d 


2 


1.21875 










3 


1.25420875420875 










4 


1.36950440777577 










5 


1.51576558282311 


0.188419 


0.468569 


-0.269131 


1.12729 


6 


1.67718058035561 


0.189303 


0.447536 


-0.213712 


1.08645 


7 


1.84742557965232 


0.194189 


0.301703 


0.217776 


0.728635 


8 


2.02343313425351 


0.197002 


0.20078 


0.545661 


0.42975 


9 


2.20346793019457 


0.197864 


0.164668 


0.672559 


0.30455 


10 


2.38615657506103 


0.198269 


0.115219 


0.715579 


0.227131 


11 


2.57169539517539 


0.198622 


0.126216 


0.821535 


0.142267 


12 


2.75870117826714 


0.198925 


0.108069 


0.897903 


0.0519525 


13 


2.9471288020509 


0.199158 


0.0926932 


0.965792 


-0.0322945 


14 


3.13672349360904 


0.199328 


0.0804934 


1.02207 


-0.10526 














oo 




0.2 










Table 5.3: Eighth moment {n — 4) of the density when F = 4. 
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N 


Mn = aN + b + c/^/N + d/N 






a 


b 


c 


d 


3 


0.829151732377539 










4 


1.13999055712937 










5 


1.45889183119874 










6 


1.78179400313294 


0.330681 


-0.23698 


-0.0197034 


0.256393 


7 


2.10668148567864 


0.326556 


-0.113869 


-0.383961 


0.55846 


8 


2.43308749295152 


0.335184 


-0.423364 


0.621536 


-0.358107 


9 


2.76069703430536 


0.335669 


-0.443708 


0.693027 


-0.428642 


10 


3.08920070504297 


0.333408 


-0.335393 


0.285749 


0.00150004 


11 


3.41837572685644 


0.33287 


-0.306358 


0.169876 


0.131419 


12 


3.74807370935471 


0.333069 


-0.318295 


0.220112 


0.0720079 














oo 




1/3 


-1/3 







Table 5.4: Fourth moment (n = 2) of the density when F = 6. 



N 


M^ = aN + b + c/^ + d/N 






a 


b 


c 


d 


3 


0.63878932696137 










4 


0.852317437834435 










5 


1.07602482170551 










6 


1.30481948267947 


0.238952 


-0.17935 


-0.0015329 


0.306508 


7 


1.5366992098612 


0.241542 


-0.256653 


0.227188 


0.116838 


8 


1.77130926929194 


0.254997 


-0.739333 


1.79534 


-1.31262 


9 


2.00817460819162 


0.253514 


-0.677215 


1.57705 


-1.09725 


10 


2.24678505124783 


0.249869 


-0.502619 


0.920542 


-0.403889 


11 


2.48676496530133 


0.249348 


-0.474543 


0.8085 


-0.278266 


12 


2.72785340315076 


0.249797 


-0.501405 


0.921547 


-0.411958 














oo 




0.25 


-0.5 







Table 5.5: Sixth moment (n = 3) of the density when F = 6. 



34 



N 


Mn = 


aN + 


b + 


c/VN + 


d/N 






a 


b 


c 


d 


3 


0.579848665870171 










4 


0.732937257370685 










5 


0.897778601877637 










6 


1.06814688756782 


0.178851 


0.00107475 


-0.227029 


0.519884 


7 


1.24220552729734 


0.18911 


-0.305123 


0.678947 


-0.231411 


8 


1.41969193536163 


0.204817 


-0.8686 


2.50959 


-1.90014 


9 


1.60002816650081 


0.201673 


-0.736901 


2.04679 


-1.44354 


10 


1.78261079422029 


0.198035 


-0.562593 


1.39137 


-0.751319 


11 


1.96700538479197 


0.198245 


-0.573912 


1.43654 


-0.801964 


12 


2.15290808717152 


0.199212 


-0.633669 


1.68802 


-1.09938 














oo 




0.2 


-2/3- -0.666667 







Table 5.6: Eighth moment {n — 4) of the density when F = 6. 



N 


Mn = aN + b + c/^/N + d/N 






a 


b 


c 


d 


3 


0.691851631655437 










4 


0.992192408359008 










5 


1.30687171332381 










6 


1.62802997072456 


0.329038 


-0.289875 


-0.448497 


0.76066 


7 


1.95072743432764 


0.302864 


0.491355 


-2.7599 


2.67751 














oo 




1/3 


-0.5 






Table 5.7: Fourth moment {n — 2) of the density when F = 8. 


N 


Mn = aN ^ b + c/Vn + d/N 






a 


b 


c 


d 


3 


0.462978101466508 










4 


0.661530611956197 










5 


0.877122153006864 










6 


1.10026340398747 


0.231106 


-0.191724 


-0.630085 


0.97549 


7 


1.32616206590665 


0.215749 


0.266649 


-1.98632 


2.10017 














oo 




0.25 


-0.75 







Table 5.8: Sixth moment (n = 3) of the density when F = 8. 



35 





Mn = aN + b + c/VN + d/N 








h 




d 


3 


n 35q3SQ45n38Q747 










4 


499919485631197 










5 


0.656534473821437 










6 


0.819712345574349 


0.165419 


0.0211561 


-0.961295 


1.19094 


7 


0.985867782818619 


0.16514 


0.0294668 


-0.985884 


1.21133 














oo 




0.2 


-1 







Table 5.9: Eighth moment (n = 4) of the density when F = 8. 

6 Conclusion 

In this work we explored two applications to the 2dOCP of the expansion of powers of the 
Vandermonde determinant based on the formalism presented in section |5J The first is the 
study of the moments of the pair correlation function on the sphere. We showed that the 
second moment satisfies an exact relation for finite number of particles A^, and we explored 
numerically the behavior of higher order moments. Also an approximation to these moments 
using a formalism based on the direct correlation function was proposed. 

The second result is on the evaluation of the distribution function of the linear statistic 
Szli l^zP" ill the soft disk geometry. The exact distribution tends to a Gaussian when A^ — )■ oo. 
We were able to compute the 0(1) correction to mean of this linear statistic for any n, and 
deduce it for a general smooth radial linear statistic. The result was checked explicitly at F = 2, 
for values of F close to 2, and numerically for F = 4 and 6. 
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Appendix 

In this appendix we present some properties of the functions 

I{k^,k2)= [[ e-'^-'H'lH'^^ dtidt2 , (A.l) 

J Jo<t2<h 

and 

J{ki, k2)= / e-'^-H'lH'^' log(max(ti, ta)) dt^dt2 , (A.2) 
Jo Jo 
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that appear in the expansion of the density around F = 2. First, it should be noticed that 



I{h,k2)+I{k2,h) = ki\k2\ . (A.3) 

Doing an integration by parts, one obtains the recurrence relation 

I{k, + 1, k^) - {k, + l)X(A;i, k^) = 2-^'~^''-\k^ + fc2 + 1)! , (A.4) 

and reiterating 

fcl! 

Similarly, for J we have 

^(fci + ^, ^2) - ^^i^ J{h^ h) = g X(fci + i, . (A.6) 

From (lA.Sp . one can obtain an alternative expression for X as a sum 

I(fa.fe) = 5:2-=-'-'(^l±^^. (A,7) 

£=0 

The asymptotic expansion of X(/ci, /C2) for large arguments, ki and /c2 of order — )■ 00, can 
be obtained by the steepest descent method. The maximum of the integrand in (lA.ll) is for 
ti = ki and ^2 = ^2- Therefore, the behavior of X(fci, ^2) will depend on whether this maximum 
is in the domain of integration < ^2 < or not, i.e., if ^2 < ki or not. If 1 ^ A;2 < ki, 
and I fci — A;2 1 / ^ 1 , then the maximum of the integrand in ( lA.lj) is deep inside the domain 
of integration, and a simple application of the steepest descent shows that X(/ci, ^2) ~ fci!/c2! • 
However for the calculations of section 1531 the behavior of X(fci, ^2) when k2 > ki is needed. In 
this situation, the maximum of the integrand is outside the domain of integration. Nevertheless, 
X(/ci, /C2) will give a significant contribution when the maximum of the integrand is "close" to 
the border, more precisely when \k2 — ki \ is of order 0{\/N). The dominant contribution to the 
integral (lA.ip will be given by the region consisting of a strip attached and parallel to the line 
ti = ^2 of width of order a/JV. To be more specific, let us do the change of integration variables 



v+ = ti + 12 and V- = ti — 12 in (lA.ip and write 



X(fci, k2) = jj^ dv^dv.e-^^^-'^-^^^^ , (A.8) 

with 

g{v+, v^) = v+-ki \og{v+ + v.) -k2 log{v+ - v^) , (A. 9) 

where V = {{ti,t2),0 < ^2 < ^i}- 9 has its maximum for v+ = v"^ = ki + k2 and f _ = f * = 
ki — ^2 (i.e. ti = ki and t2 = /c2)- Expanding g to the second order around its maximum we 
obtain 

Iik„k2) ~ e-^[(--<)^-(-^^)1-«(--^^)(-^^)ci.,rf._. (A.IO) 



37 



Although the domain of integration V cannot be simply expressed in terms of the variables u+ 
and f_, due to fast Gaussian decay of the integrand in flA.lOp . the dominant contribution is 
indeed obtained by the strip along the line ti = t2 mentioned above. With this in mind, it is 
clear that one can extend the domain of integration for t>+ to ]— oo, +00 [ and for t>_ to [0, +oo[, 
up to exponentially small corrections. Then, performing first the integral over f+, we obtain 



Xik^.k^) ^ e-^'-^^k\'kl\ / e" 2(^1+^2) dv_ . (A.ll) 

y ki + k2 Jo 

Performing now the integral over f _ gives 

I{ki, k2) ~ e-^^-'=2A;t^A;2^^7rv/M^erfc ( ~ | , (A.12) 

\^2{ki + k2) ) 

where erfc(x) = 1 — (2/ a/tt) e~*^ dt^ is the complementary error function. Recalling Stirling's 
formula for the factorial, (IA.12P can be written as 

- ierfcf^i^^l . (A.13) 



ki\k2\ 2 \^j2{ki + k2) 
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